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Abstract 

Earlier attempts to extract parameters of kinetic freeze-out in Pb+Pb collisions at 158AGeV 
are critically reviewed. Many of these analyses have used approximations which have sig- 
nificant impact on the extracted parameters. Simple estimates are obtained which attempt to 
avoid the most critical approximations. It is pointed out that constraints based on pion inter- 
ferometry are less reliable than those from momentum spectra. A universal set of freeze-out 
parameters from transverse mass spectra would require T > 135 MeV and (fir) < 0.35. 



1 Introduction 



Experiments with ultrarelativistic heavy ions at the CERN SPS which are on the 
quest for a new state of strongly interacting matter have obtained a wealth of in- 
formation during the last years (see e.g. [1,2]). While there are several interesting 
non-trivial observations in these experiments, our understanding of the evolution of 
the reactions is far from complete. However, for an interpretation of the promising 
possible signatures of the quark-gluon plasma a better understanding of the reaction 
dynamics is mandatory. 

Hydrodynamical models or simple parameterizations have been used as tools to 
help understanding some of the issues in heavy ion reactions. While even the more 
complex versions of such models are believed to be crude approximations, it is 
worthwhile to use them and also the simpler parameterizations to test strongly 
different scenarios and provide first guidelines on the possible dynamics, if the 
reaction mechanisms are not too far away from the concept of local thermal equi- 
librium. One of the key questions in this context is the partition of the excitation 
energy into collective and thermal (i.e. random) motion of the produced particles. 
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This is mostly characterized by the average values of the temperature and the ra- 
dial flow velocity of the finally observed hadrons, the so-called kinetic freeze-out 
parameters. 

Still, in such comparisons one has to be cautious to 

(1) use only those experimental results which are most likely too share the same 
equilibrium conditions and 

(2) avoid approximations that are wrong even in the most optimistic cases. 

In the following I will study some of the experimental results on momentum spec- 
tra and interferometry that have been used in applications of these hydrodynamical 
models and check their implications for the estimates of freeze-out parameters. I 
will first use estimates based on effective slopes of transverse momentum spectra 
for different particle species and compare them with the results of hydrodynamical 
parameterizations of neutral pion spectra. I will then turn to the interpretation of 
pion interferometry results. There the momentum dependence of effective radii ex- 
tracted from fits to two-particle correlation functions is used to extract information 
about temperature and flow. As the method to obtain this information is much more 
indirect, it is naturally more complicated to interpret the corresponding results. I 
will try to demonstrate that some of the recipes established should be critically 
revisited. 



2 Analysis 

2.1 Inverse Slopes 

The transverse momentum spectra of hadrons are believed to be frozen in a rela- 
tively late stage of the nuclear collision. If there is local equilibration, the elemen- 
tary fireball volumes can be associated with local temperatures. In a finite system 
there will necessarily be hydrodynamic motion because of the pressure gradient. 
Therefore, at freeze-out the system may be characterized by local temperatures and 
transverse flow velocities which will determine the shape of the transverse mo- 
mentum spectra.0 If the composition of the system is relatively homogeneous, the 
freeze-out may occur at a relatively similar temperature throughout the whole vol- 
ume. It is therefore reasonable to speak of one general freeze-out temperature. One 
should keep in mind, however, that this is a crude approximation: It is more a ques- 
tion of how different all the local temperatures are and for which particle species 
they are relevant. This simplification may be used for the discussion of momentum 

2 There is an additional dependence of the shape on particle production by resonance decay 
which I will not discuss in detail here. 
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spectra, where a superposition of distributions with different temperatures may well 
be approximated by one distributions with an average temperature. This simplifi- 
cation is much more crucial in the case of pion interferometry, as I will discuss 
below. 

The general idea in using different particle species to estimate these parameters 
relies on the fact that in a system with a temperature and a collective flow veloc- 
ity, heavier particles will acquire higher momenta from the shared velocity because 
of their larger mass and will exhibit flatter momentum spectra, i.e. larger inverse 
slopes. In [3] the authors have presented inverse slopes T cS of transverse mass spec- 
tra of different particle species and used those to estimate a freeze-out temperature 
and flow velocity from the relation: 

T cS = T + m-(p T } 2 . (1) 



They obtain a freeze-out temperature of T pa 140 MeV. This analysis has been 
criticised as neglecting relativity [4,5], the authors of the latter publications present 
a better estimate of the slopes: 



' . ii = v ' - (Pt) 2 T = ■ (2) 
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Then the limit of zero particle mass: 
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is used to reestimate the freeze-out temperature. They obtain T pa 90 MeV using 
the data of [3]. Both groups use also somewhat more sophisticated models to obtain 
freeze-out parameters which support their simple temperature estimates. 

I would like to demonstrate the ambiguities in such estimates by using the cor- 
rect relativistic formula 2 to calculate freeze-out temperatures from the same slope 
parameters. For this purpose I have to assume pr = (jPr)- This is of course an ap- 
proximation which should however not be as severe as those used in [3] and [4,5]. 
The slopes and the values of (jp T ) used are given in table 1. It turns out that not all 
three particle species investigated can be described with the same set of parame- 
ters. Comparing positive kaons and protons yields T = 139 MeV and (f3 T ) = 0.33, 
while for negative kaons and antiprotons one gets T = 152 MeV and (f3 T ) = 0.28. 
However, estimates using pions and (anti)protons yield T pa 100 MeV and {(3 T ) ~ 
0.41 and those for pions and kaons yield T pa 84 MeV and (f3 T ) ~ 0.51. 

Although the authors of [3] exclude the region of m T — m < 0.25 GeV/c from 
the fit for pions, it is still very likely that the slopes for pions are influenced by 
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Fig. 1. Relations between temperature and flow velocity extracted from hydrodynamical 
parameterizations of neutral pion spectra from WA98 (grey band) and from slope parame- 
ters of proton and antiproton (NA44 and NA49, solid lines), kaon (NA44, dashed lines) and 
deuteron spectra (NA44 and NA49, dotted lines). The light grey area shows the allowed 
region from pion interferometry analysis by NA49. (For details see text.) 

resonance decays which would result in a lowering of the inverse slope and thus 
implicitly lower the extracted freeze-out temperature. Kaons and protons should be 
relatively free from such influence, so that the estimate using these particle species 
should be more reliable. The temperature T = 140 — 150 MeV is also in line with 
the hydrodynamical fit analysis in [3]. 

The detailed analysis in [4] which obtains lower freeze-out temperatures uses addi- 
tional information, namely the normalization of the spectra and two-particle corre- 
lations. It is known from other analyses (see e.g. [6]) that two-particle correlation 
which can provide constraints on (Pt) 2 /T within hydrodynamical models effec- 
tively gives very little information on the temperature but constrains mostly the 
flow velocity (see light grey area in figure 1). On the other hand the normalization 
of spectra of different particle species is related to chemical freeze-out, while the 
spectral shape is influenced by kinetic freeze-out ^\ To assume that these two cor- 
responding temperatures are the same, as was apparently done in [4], is very likely 
wrong.0 

From the large pion-nucleon cross section one would assume that these two parti- 

3 There is an influence of the chemical temperature on the spectral shape via the resonance 
decays - this should however be relevant only for pions. 

4 A similar assumption is used in the analysis in [7,8] which will not be discussed here. 
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cle species are most likely to share the same freeze-out conditions. However, the 
influence of resonance decays in the pion spectrum forbids the simple use of slope 
parameters to estimate freeze-out conditions. To perform such a comparison the 
slope parameters for heavier particles will be used together with a full hydrody- 
namical parameterization of neutral pion spectra as described in [9]. The hydrody- 
namical computer program described in [10] has been modified to describe neutral 
pions and accommodate a Woods-Saxon density profile in addition to the default 
Gaussian profile and was used to obtain best parameter sets for T and (Pt)- The 
2<7-ailowed values obtained for different profiles (Gaussian and Woods-Saxon of 
different diffuseness) are shown in figure 1 as a grey band. Among these fits, those 
obtained with large diffuseness (as e.g. the Gaussian) tend to yield smaller flow 
velocities for the same temperature but also constrain the allowed parameters to the 
high T region [9]. In addition, the slope parameters for kaons, protons, antiprotons 
and deuterons from NA44 and NA49 as summarized in table 1 have been used with 
the help of equation 2 to obtain similar relations between T and (Pt)- All these 
estimates are also displayed in figure 1. 

The NA49 collaboration has also published freeze-out constraints obtained from 
negative hadron and deuteron transverse mass spectra [6]. The negative hadrons 
are dominated by pions which should allow to use them in hydrodynamical descrip- 
tions assuming the pion mass, however, the analysis was performed using a non- 
relativistic approximation and ignoring resonance production. The deuteron spectra 
are not affected by resonances but still suffer from the non-relativistic approxima- 
tion. The corresponding parameters are therefore not reliable. In fact, the estimate 
from the slope included in figure 1 differs considerably from the constraints in [6], 
esp. at large velocities. While the curves for kaons, protons and deuterons seem to 
converge around T = 140 MeV and (Pt) = 0.33, the neutral pion results have a 
further bias to somewhat higher temperatures. One should note that the analysis of 
the pion spectra differs from the others, as 

(1) a full parameterization including the effects of resonance decays has been 
performed compared to using an effective slope together with an approximate 
relation connecting T and (p T ) and 

(2) they cover a different m T range - the high m T tail of the neutral pion spectrum 
has also been used as an upper limit as described in [9]. 

Still this figure suggests that a common parameter set for all the momentum spectra 
considered here can only be obtained for T > 135 MeV and (Pt) < 0.35. This 
appears to contradict constraints obtained from pion interferometry measurements 
[6] (see light grey area in figure 1), which favor larger transverse flow velocities. 
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2.2 Interferometry 



It is well known that correlations between position and momentum space modify 
the radii measured in pion interferometry. One of the most interesting examples of 
such a correlation is collective flow [14,15], both longitudinal and transverse, al- 
though there are other examples, like a source with local temperatures T = T(x^), 
as I will discuss below. In the idealized case of pure collective motion (no thermal 
contribution) only the source element moving towards the detector will be mea- 
sured, so the radius extracted from interferometry will correspond to this small 
elementary volume. If there is additional thermal motion, the probability of a parti- 
cle to still be measured while being radiated from another source element at another 
space point and directed slightly away from the detector will increase. The effective 
source size will therefore decrease, if the amount of collective motion relative to 
thermal motion is large, and vice versa. In addition, this effect is a function of the 
momentum of the particle, and in most scenarios the effective radii will decrease 
with increasing transverse mass of the particle. 

Most quantitative studies of these effects start from the Cooper-Frye formula [16] 
for the momentum distribution of the source: 

d 3 N _ 2J+1 r 3 f_ 

dp 3 (2vr) 3 J ^ exp { [p^ — ji] /T} =F 1 ' ( ' 



In this formula the factor p^ in the numerator and the exponent p^v M — fi originate 
from the Lorentz transformation with v M = v being the local four- velocity of 
the source elements. In general also T is again a local variable. // is the chemical 
potential and the =f signs are relevant for bosons and fermions, resp. 

There has been substantial work on extracting expressions for observables in in- 
terferometry, the most systematic approach using a parameterization of the source 
distribution and a saddle point approximation to extract approximate analytic re- 
lations for the effective radii [17]. A relatively general result is obtained using a 
source distribution expressed as: 

K-u x 2 + y 2 (r — r ) 2 rf 
~T 2~R| 2(At) 2 ~ 2(A V )\ 



S(x,K) = 



m t ch(r] — Y) 



cxp 



Gaussians are used for the space-time and rapidity distributions, the momentum 
spectrum is approximated as a Boltzmann distribution: 

d 3 N ( K-u\ 
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For a model with transverse expansion, the authors obtain approximate relations 
for the transverse and longitudinal radii: 



Rrp ~ Rq 
Rl ~ T~0 



-1/2 



^ ch( ^- F) -d?(^) + (A^ 



-1/2 



(7) 
(8) 



For a boost invariant longitudinally expanding source the expressions are more 
simple, e.g. for the longitudinal radius one obtains: 



T 



Rl ~ T J . 

V VtlT 



(9) 



As it turns out [6] that the source rapidity and the pair rapidity are similar i] m Y, 
also the expression for the transverse radius may be simplified: 



Ri 



Rg 



1 + — • m T 



1-1/2 



(10) 



These two latter equations are mostly used to extract information from experimental 
two-particle correlation functions. There are a number of assumptions which have 
to be made to obtain equations 7 and 8 and even more for equations 9 and 10. Some 
of these are discussed in [17], however there are further assumptions which I will 
try to concentrate on here. 

To be able to use equation 10 to extract as in [6], one of the major assump- 

tions is that the parameter R G , which contains information on the geometrical size 
of the source, does not depend on the transverse momentum of the emitted parti- 
cles. This is not necessarily correct. In general, there will not be a universal value 
for the temperature, it will be a local parameter depending on the position of the 
source element with respect to the center or the surface as discussed above. 

This induces an additional correlation between momentum and coordinate space. If 
the variations in T are very strong, this would of course invalidate all the attempts to 
extract the freeze-out temperature. If the variations are relatively weak, it may still 
be reasonable to use one effective temperature. It may e.g. be used to characterize 
momentum spectra, as for those the spatial distribution and therefore also these 
correlations are not directly relevant. 

However, they are crucial to the analysis of pion interferometry, especially to the 
above mentioned method of extracting Consider e.g. that a certain fraction 

of pions is emitted before freeze-out, from the surface of the source, and that this 
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Fig. 2. Effective source size as a function of the transverse mass of the particle for a super- 
position of two different sources (see text). 

fraction increases with increasing momentum. One might try to describe this with a 
superposition of two sources. Just as one example one may consider the following 
- arbitrary - scenario: One source of effective size (Gaussian) i? x = 10 fm and 
temperature T\ = 200MeV and a second with R 2 = 2 fm and T 2 = A00MeV.^\ 
If one investigates produced particles in different intervals of transverse mass, the 
relative contributions of these two sources will change, and implicitly the effective 
source size as seen in interferometry. This is shown in figure 2, where the effective 
source size - again approximated by a Gaussian - is shown as a function of the 
transverse mass for the case where the two sources discussed above both contribute 
the same fraction to the total multiplicity of particles. It can be seen easily that 
the effective source size decreases significantly with increasing transverse mass. 
Of course, this scenario has not been derived from physical arguments here, so it 
should only serve as an illustration of possible effects of such correlations between 
momentum and coordinate space. How large such effects might be would require 
a detailed model calculation which is not the scope of this paper. One should only 
note that in general Rq in equation 7 is expected to be a function of mx within the 
model used. Values of ft\jT should therefore be considered as upper limits unless 
more detailed model comparisons are performed. Pion interferometry is much more 
sensitive to such effects than momentum spectra. 

A more fundamental caveat to the use of the tut dependence of HBT to estimate 
freeze-out parameters relates to another approximation used in extracting the ana- 
lytical formula which are applied in most of the analyses. Specifically, the authors 
of [17] use the Boltzmann approximation to the distribution function. While this 
may be a good approximation for most of the pions used in experimental HBT 
analyses, it is well known that the deviations to the true Bose-Einstein distributions 



For simplicity exponential distributions in transverse mass are assumed. 
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Fig. 3. Inverse slopes (equation 11) relative to pure exponential slopes as a function of 
energy for different distribution functions. 

are crucial at low momenta. This is illustrated in figure 3, where the inverse slopes 
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eff 



dE 



log f(E) 



of Bose-Einstein and Fermi-Dirac distributions 



E d 3 N ^ 1 

dp 3 exp(-f)^l' 



are shown. Strong deviations from Boltzmann behaviour are seen for both distribu- 
tions for E /T < 1. 

As in the analyses in question most of the m T dependence is seen at relatively low 
momenta, one should carefully reinvestigate the interpretation of these data with 
the correct statistical distributions. The trend of such a corrected analysis may be 
guessed from the behaviour of the slopes. Using the crude assumption that the BE- 
distributions may be locally approximated by a Boltzmann distribution, one might 
then use the formalism developed in [17] to estimate the m T dependence. However, 
now the effective slope for a Bose-Einstein distribution (T BE ) as shown in figure 3 
should be used instead of the Boltzmann temperature T. In figure 4 the transverse 
radii as a function of transverse mass are shown for three different estimates: 

(1) Following the formalism as in [6], i.e. using the Boltzmann-like temperature in 
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Fig. 4. Transverse radii for pion interferometry as a function of transverse mass according to 
equation 10. The dashed line uses the Boltzmann temperature with values of T = 120 MeV 
and (fir) < 0.55, following the results of [6]. The solid line uses the same parameters with 
a correction for the effective slope of a Bose-Einstein distribution. The same correction is 
used for the dotted line but T = 85 MeV is assumed. 

equation 10 and the results obtained in the article, i.e. T = 120 MeV, {(3r) = 
0.55 and Rq = 6.5 fm, 

(2) substituting the temperature by the corresponding effective slope of a Bose- 
Einstein distribution, but using the same parameters as above and 

(3) substituting the temperature by the corresponding effective slope of a Bose- 
Einstein distribution, but now using T = 85 MeV, (fir) = 0.55 and Rq = 
7.1 fm. 

It can be seen from this example that cases 1 and 2 show significant differences, 
and that for a rough agreement with 1 the temperature in the corrected version has 
to be modified from 120 MeV to 85 MeV as in case 3. 

In this context it is interesting to investigate also results for proton source sizes. 
While the mechanism generating the correlation of particle pairs is different com- 
pared to pions - a strong interaction quasi-resonance instead of Bose-Einstein statis- 
tics, the picture of effective source sizes modified by the interplay of random ther- 
mal and collective hydrodynamical motion remains valid. People have therefore ap- 
plied similar ideas to pp-correlation radii. The experiments NA49 [19] and NA44 
[20] have measured proton correlations at midrapidity in Pb+Pb collisions - they 
obtain Gaussian radii in the range of R pp = 3.5 — 4.0 fm which are considerably 
smaller that those extracted from pion interferometry (R n7T = 6 — 7 fin). Similar 
values have been estimated by applying coalescence models to measured proton 
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and deuteron multiplicities [20] . The smaller effective radii have been attributed to 
the much larger transverse mass of the protons, if a scaling with 1 / ^Jmr which may 
be derived as an approximation from equations 9 and 10 is used for both protons 
and pions. 

From the above argument on BE-distributions it is obvious that equivalent modi- 
fications are important here. However, the effective slope for a FD-distribution is 
larger at low momenta. In addition, the exponent of the FD-distribution contains a 
non-vanishing chemical potential which makes it much more likely that the expo- 
nent is small and that modifications are important. 

Actually, another important argument is relevant here. If, as commonly believed, 
chemical freeze-out occurs relatively early before kinetic freeze-out, the relevant 
picture for the latter is the canonical ensemble, as the particle number for each 
species is fixed, so that there is no chemical potential as a parameter. The relevant 
energy in the statistical distributions is the kinetic energy and not the total energy, 
as no new particles are produced and the ground state energy of each particle is the 
rest mass. The Cooper- Frye formula should then be changed to: 



The major change for the transverse degrees of freedom would be to use mr — rn 
instead of m T . This should hold for both bosons and fermions, but while it might 
have only small effects on the pions, it will be dramatic for protons because of 
their large mass. It is beyond the scope of this paper to estimate the effect of such 
modifications on the m T spectra in general or, more specific, on the m T dependence 
of proton source sizes. However, it is also premature to conclude that proton and 
pion source sizes are compatible without performing detailed calculations. As the 
values of m T — m for protons and pions are roughly similar in the experiment, 
no difference in the radii would be expected between protons and pions from their 
transverse masses. This would however contradict the measured source sizes. 



3 Conclusion 

I have tried to illustrate above that the freeze-out parameters in ultrarelativistic 
heavy ion collision are not well established quantitatively. One should be careful 
to use only data sets which very likely share the same conditions. Also, the used 
approximations should be thoroughly investigated. The simple re-analysis of some 
of the available data performed here indicates the following: 

• Transverse momentum spectra of deuterons, protons, kaons and pions seem to 




(13) 



11 



indicate that compatible parameter sets can be obtained only for T > 135 MeV 
and {fa) < 0.35. 

• The m T dependence of interferometry radii has to be reinvestigated using true 
statistical distributions instead of Boltzmann approximations at low momenta. 
This would also require to study the role of modifications related to the change 
at chemical freeze-out to a system with fixed particle number. 

• Other possible sources of correlations between momentum and coordinate space 
may cause further alterations of the tut dependence of interferometry radii. 

• The difference in pion and proton source sizes is very likely not explained by the 
difference in transverse mass of the two particles. 

The three last observations shed doubt on the use of pion interferometry as the main 
ingredient to obtain constraints on the amount of transverse flow at freeze-out. A 
more thorough analysis of the pion interferometry results is necessary to clarify the 
situation. 

Two possible scenarios may be envisaged which explain the data discussed here: 

(1) Pion HBT does not yield reliable constraints on kinetic (thermal) freeze-out.Q 
It may actually be that pion interferometry measures a different, later time than 
proton correlations and momentum spectra. This late "wave function freeze- 
out" may be related e.g. to very soft or mean field interactions of pions which 
do not alter the momentum distributions as mentioned in [21]. Momentum 
spectra, however, favour T > 135 MeV and ((3 T ) < 0.35. 

(2) Freeze-out does not occur with one even approximately universal temperature 
for different species and also different parts of the source. 
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Table 1 

Parameters from transverse mass spectra of pions, kaons and protons in central 158AGeV 
Pb+Pb collisions. The slope parameters were taken from [3,13,11,12], the average trans- 
verse momentum has been calculated assuming an exponential shape of the m-f spectra in 
the fit range used. 
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